Transport control by coherent zonal flows in the core/edge transitional regime 
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3D Braginskii turbulence simulations show that the energy 
flux in the core/edge transition region of a tokamak is strongly 
modulated - locally and on average - by radially propagat- 
ing, nearly coherent sinusoidal or solitary zonal flows. The 
flows are geodesic acoustic modes (GAM), which are primar- 
ily driven by the Stringer- Winsor term. The flow amplitude 
together with the average anomalous transport sensitively de- 
pend on the GAM frequency and on the magnetic curvature 
acting on the flows, which could be influenced in a real toka- 
mak, e.g., by shaping the plasma cross section. The local 
modulation of the turbulence by the flows and the excitation 
of the flows are due to wave-kinetic effects, which have been 
studied for the first time in a turbulence simulation. 

a. Introduction — It is now commonly believed, that 
the transport in the tokamak core is controlled by zonal 
flows 0-|§]. In the plasma edge, the flows were not stud- 
ied thoroughly yet, but they tend to be weak 0] (al- 
though they can possibly completely quench the tur- 
bulence pI). The zonal flows in the core and at the 
edge were found to be incoherent, random fluctuations 
[p|,p|j^. In contrast to this, in the transitional regime be- 
tween core and edge, strong radially coherent sinusoidal 
or solitary zonal flows are ubiquitous according to the 
numerical simulations described below. The zonal flows 
in the transitional regime are essentially geodesic acous- 
tic modes M, i.e., the poloidal rotation is coupled to an 
{m,n) — (1,0) pressure perturbation by the inhomoge- 
neous magnetic field, which results in a restoring force 
and hence an oscillation. Transcending the present mod- 
els of the zonal flow generation based purely on Reynolds 
stress ||l[, the major part of the flow energy in the tran- 
sitional regime is apparently generated by the Stringer- 
Winsor term p|, i.e., the torque on the plasma column 
caused by the interaction of pressure inhomogeneities 
with the inhomogeneous magnetic field. An order of mag- 
nitude variation of the shear flow level and average trans- 
port has been observed upon a sole modiflcation of the 
curvature terms acting on the flows, at flxed curvature 
terms acting on the turbulence. The pressure inhomo- 
geneities are driven by anomalous transport modulations, 
which can be understood by a drift-wave model for the 
wave- kinetic turbulence response to the flows. Some pre- 
dictions of the model have been verified subsequently in 
numerical experiments. 

Zonal flows in the transitional regime are GAMs — 
The numerical turbulence simulations have been carried 
out using the three dimensional electrostatic drift Bra- 
ginskii equations with isothermal electrons, including the 
ion temperature fluctuations with the associated polar- 
ization drift effects, the resistive (non-adiabatic) parallel 
electron response and the parallel sound waves (a subset 



of the equations of Ref. ||]). The nondimensional param- 
eters have been varied around a reference parameter set 
resembling the transitional core/edge regime: ad — 0.6, 
e„ = 0.08, q — 3.1, r = 1, 77^ = 3, s = 1. The compu- 
tational domain is a flux tube winding around the toka- 
mak for three poloidal connection lengths. The radial 
and poloidal domain width is SOLrb, with the resistive 
ballooning scale length, Lrb ■ For a definition of these pa- 
rameters and units see Ref. H . The parameters are con- 
sistent with the physical parameters i? = 3 m, a = 1.5 m, 
i„ = 12 cm, n = 3.5 X 10^^ m^^^ z^s = 4, Bq = 3.5 T, 
and T = 200 eV, Lrb = 3.6 mm, ps = 0.82 mm. At 
these parameters the ITG mode turbulence is the domi- 
nant cause of the heat flux M , and the restriction of the 
domain to a flux tube is justifled |0]. 

Viewed as a function of radius and time, the flux sur- 
face averaged poloidal E x B flows [fig. |l| (a)] start as 
an irregular pattern of radially propagating independent 
wavelets and merge into a radially coherent standing 
wave later on. The final standing wave pattern consists 
of GAM oscillations. They can be described by suitable 
poloidal Fourier components of the vorticity and pressure 
evolution equations (neglecting the parallel sound wave) , 
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with (.) denoting the fiux surface average, the pressure 
fiuctuations p = n + j^—Ti, the poloidal fiow velocity 
ve = dx4>j ^iid the source terms s^/sp of fiow/pressure 
due to Reynolds stress/anomalous transport, respec- 
tively. For easier reference, the two different curvature 
terms in the equations have been adorned with the factors 
Ci and C2 , which are both 1 in the turbulence equations 
for a low aspect ratio circular tokamak. The curvature 
term Ci is the Stringer- Winsor term, the term C2 repre- 
sents the up-down asymmetric compression of the plasma 
due to the poloidal rotation. For all parameters used in 
the numerical simulations, the zonal flows oscillate in the 
stationary state with a frequency within 5% of the eigen- 
frequency of Eqs. (pJj3) without source terms and with 



the approximation 
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Cs/R. Since the 



10t)/(3 + 3t)cJR 
parallel sound frequency Cs/(qR) is much lower than the 
GAM frequency, the neglect of the parallel sound wave 



is justified. The energy balance equation of the GAM 
oscillations is according to MF 
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The average contributions to the GAM energy from the 
Reynolds stress term (s^) and the Stringer- Winsor term 
{Sp sin 9) have been listed in table | for varying turbulence 
parameters. In the transitional regime with its strong co- 
herent zonal flows most of the flow energy is generated by 
the Stringer- Winsor term, while with decreasing temper- 
ature toward the very edge the Stringer- Winsor energy 
input eventually becomes negative indicating a braking 
force on the flows, while simultaneously we get weak in- 
coherent flows. One is tempted to attribute the decrease 
of the zonal flows towards the very edge to the Stringer- 
Winsor term. Indeed, eliminating the flow source term 
(sp sin 0) due to the anomalous transport in the numer- 
ical simulations leads to relatively strong coherent flows 
even for parameters in the very edge. 

With the natural drive of the GAM being apparently 
the Stringer- Winsor term, it has been found that altering 
the amplitude of the curvature terms Ci , C2 acting on the 
flows but keeping the curvature terms acting on the tur- 
bulence modes fixed, changes the flow amplitudes and the 
anomalous transport by one order of magnitude. Empir- 
ically, the flow level rises with increasing Ci and decreas- 
ing uj. One reason is, that the higher Ci/w ~ A/C1/C2 
is, the higher is the contribution of the anomalous trans- 
port source term Sp to the flow energy (^. At constant 
ratio Ci/uj the flows are still somewhat increasing for 
decreasing uj. This is understood, since with increasing 
oscillation period the flows have more time to influence 
the turbulence resulting in an increased Sp. 

The cause of the pressure perturbations driving the 
flows are local modulations of the anomalous transport 
(the usual Stringer spin-up mechanism is ineffective due 
to the relatively long sound transit time in the con- 
sidered regime). Plots of the radial pressure transport 
{vrp) — Q{t, t) and its up-down antisymmetric compo- 
nent {vrPSvuO) — U{r,t) are shown in fig. |] (b) and (c). 
(Note that the pressure source term is the divergence of 
the anomalous pressure flux, i.e., Sp = —dr{vrp), and 
(spSin^) — —drU.) In the initial phase of flow genera- 
tion, U develops dipoles around the flows [flg. |] (d)], in 
which U has always the same sign as the local shear- 
ing rate. These dipolar transport structures generate 
the pressure up-down asymmetries driving the GAMs. 
As soon as sufficiently strong flows exist, the anoma- 
lous transport Q develops a striking peaking at the radii 
of positive flow resembling transport fronts propagating 
with the flow "waves". In addition to the dipole struc- 
tures U develops a unipolar part component, whose sign 



depends on the propagation direction of the correspond- 
ing flow. These up-down antisymmetric transport fronts 
can be viewed as avalanches running outward on the 
lower half of the torus {9 < 0) and inward on the up- 
per half. 

The unipolar up-down asymmetries have been found 
to be responsible for the setup of the flow pattern. If in a 
numerical experiment the GAMs are initially set to zero 
outside of one flow peak, the turbulence is still capable of 
moving this flow into the original direction, until a new 
standing wave pattern has been formed. The necessary 
radial GAM energy flow due to the turbulence has been 
found to be primarily caused by the unipolar up-down 
asymmetries. Similar events happen, if a numerical sim- 
ulation is started with a GAM pattern with the wrong 
wave-length from a simulation run with different GAM 
parameters. In that case, strong unipolar up-down trans- 
port asymmetries develop, which enforce the equilibrium 
flow propagation velocity, i.e., flow wavelength. 

The described peculiar modulation of the transport by 
the flows is absent for weak diamagnetic drift and vanish- 
ing gyro radius, such as in the resistive ballooning regime 
[Q. Instead, the shear flows simply weaken the turbu- 
lence. Apart from this, there is a tendency to flatten 
pressure gradients, i.e., to eliminate the pressure fluctua- 
tions associated with the GAM, resulting in the observed 
braking force. 

The wave-kinetic effects The dependence of the 
transport modulation on drift effects and the gyro ra- 
dius suggests a simplified drift wave model containing 
only the radial mode coupling due to the polarization 
drift, eliminating two fluid and curvature effects. Since 
in the numerical studies the mode wavelengths are not 
small compared to the zonal flow scales, we refrain from 
a geometrical optics approach m and instead move back 
to the linearized adiabatic drift wave equation. 
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with Dt — dt + VE{x,t)dy. We study the impact of the 
polarization drift term —Dtp^A±n up to flrst order in p^, 
without assumptions on the ratio of flow vs. turbulence 
scales. The 0-th order time evolution (with y in Fourier 
space) is 



4>oix,t) = (t>i{x)'4){x,t) 
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with initial amplitude 4'i{x) and the flow-and-drift in- 
duced phase factor 



^p{x,t) = exp[-iky^], £, ^ [ve{x,t) + ad] 

Jo 
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Inserting (]m into (h) results in the flrst order correction 
(/)i due to p^ 

4>i{x,t)ijj*{x,t) = -ikyttdpl I i^*{x,T)A±{4)itl;{x,T))dT+ 
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in which [f{. . . , t)]1% = /(■ ..,<)-/(..., 0). From this, 
the change in turbulence intensity can be computed to 
first order, 
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with a suitable local kx{x,T) = ^kyS,' {x^t) + Im((/)^/0j). 
The term involving the diamagnetic drift velocity a^ cor- 
responds to the advection of mode intensity in response 
to the shearing distortion, while the other term is analo- 
gous to the adiabatic compression of a wave field. 

If the initial (j)i has no radial structure {k^ = 0) and 
the shear flows do not change with time, we obtain due 
to radial mode advection 
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which explains the empirical peaking in turbulence inten- 
sity at the locations of positive flows. To verify the pres- 
ence of the advection term, a complementary scenario has 
been simulated, with the turbulence initially conflned to 
a small region and a constant linear shear flow ve o(. x 
enforced upon it. This results in a motion of the turbu- 
lence maxima towards increasing ve , which confirms the 
presence of the advection term (see fig. 0). In another nu- 
merical experiment, a stationary bell shaped flow profile 
has been superposed on an initially radially homogeneous 
turbulence field, which has been prepared with the zonal 
flows switched off. After the turbulence rises transiently 
at the flow maximum, it drops to a level below the initial 
one. This corroborates that the turbulence amplification 
at the flow maxima is not due to an increase of drive at 
the flow maxima [|lO| , but due to a transient wave-kinetic 
concentration of fluctuation energy. 

The numerically observed radial dipole layers of up- 
down antisymmetric transport apparently result from the 
compressional terms in (^ acting on the mode structure 
enforced by the magnetic shear, namely 
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Hence, with s = 1 a shear flow with positive v'^ reduces 
fc^ on the upper side of the tokamak {0 > 0) increasing 
the mode amplitude there, while the mode is attenuated 
for (9 < 0). The dependence of the shear flow action 
on the magnetic shear has been verified in a numerical 
experiment with s = —1, in which the modes are am- 
plified for 0Dg < 0. Moreover, a "swinging through" 
effect hast been observed, i.e., the turbulence intensity 
decreases again, when fc^ rises after k^ has gone through 
zero. 

Last, we consider the unipolar up-down transport 
asymmetries induced by a propagating positive shear flow 



(which is accompanied by a general transport peaking 
due to the advective effect). With VE[x,t) = V{x — vt) 
(propagation velocity v) we obtain from definition (M) 
\'[x,t) = -v-'^V{x-ut), i.e. 

fc,(x, tf = (fc,o " kyS,'? = kl{se + u-^V{x - ut)f. 

(12) 

A reduction of fc^ and corresponding amplification of the 
turbulence modes via eq. (||) occurs for sOvV < 0. The 
sign of this effect of a moving flow has been confirmed 
numerically for positive and negative shear. 

Conclusions — GAMs, oscillating zonal flows, have 
been found to be the main mechanism controlling the tur- 
bulence level in the transitional core/edge regime. Quasi- 
stationary zonal flows are unimportant in the simula- 
tions, because of the strong restoring force due to the 
pressure imbalance generated by the magnetic field in- 
homogeneity, which cannot be short-circuited along the 
magnetic field lines due to the long parallel sound transit 
time in the transitional regime. Primarily, the flows are 
not driven by Reynolds stress but by the pressure asym- 
metries on a flux surface generated by modulations of 
the anomalous transport. These modulations in turn are 
caused by the flows, on one hand by the radial advec- 
tion of turbulence energy concentrating the turbulence 
in locations of increased flow in electron diamagnetic di- 
rection, and on the other hand by adiabatic compression 
effects on the wave field which together with the mag- 
netic shear result in up-down asymmetries of the trans- 
port. The peculiar nature of the drive mechanism leads to 
propagating peaked flow structures acompanied by pro- 
nounced transport fronts, which sometimes come close to 
propagating solitons. We have studied for the first time 
the action of wave-kinetic effects in a numerical simu- 
lation. Similar results should be expected for the core, 
except that the polarization due to the ion larmor ra- 
dius has to be replaced by the neoclassical polarization 
due to the banana width [[ll| . I.e., we expect the wave- 
kinetic transport modulation in the core to be signifi- 
cantly stronger than in in the edge. GAMs on the other 
hand are less important in the core due to the shorter 
parallel connection length {q ^ 1). 

On one hand, the GAM drive efficency depends on the 
nature of the turbulence, in that it depends on the pres- 
ence of finite ps effects. Absence of these, such as in 
the resistive ballooning regime leads to a strong damp- 
ing of the GAM due to the anomalous diffusion erod- 
ing the pressure fluctuations connected with the GAM. 
However, apart from the turbulence the zonal flow am- 
plitude is influenced strongly by the linear properties of 
the GAM itself, such as its frequency or the torque ex- 
certed on a pressure fluctuation. Manipulation of both of 
them can lead to a reduction of the transport by up to an 
order of magnitude. Consequentially, any discussion of 
the anomalous transport focusing on the influence of the 
magnetic geometry on the turbulence drive falls short of 
an essential factor, if the influence of the magnetic geom- 



etry on the flows is neglected. Since the GAM properties 
could be influenced in a real tokaniak by, e.g., shaping the 
plasma column, the coherent flows should be taken into 
consideration to further reduce the transport in advanced 
tokamaks. Moreover, the clear signature of the radially 
coherent flows in the numerical simulations makes them 
an interesting target of experimental investigation, e.g., 
by means of microwave reflectometry. 

This work has been performed under the auspices 
of the Center for Interdisciplinary Plasma Science, a 
joint initiative by the Max-Planck-Institutes for Plasma 
Physics and for Extraterrestrial Physics. 
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FIG. 1. Time evolution of (a) poloidal E x B flow pro- 
file, (b) pressure flux profile (vrp), (c) up-down asymmetric 
pressure flux {vrP sin d) for the standard parameters; (d) in- 
stantaneous profiles at i = 80 of the flow (above), pressure 
flux (solid, below), and the up-down asymmetric pressure flux 
(dashed) . Note the dipole layers of up-down asymmetric pres- 
sure flux around the flows, e.g., at x = —7. The unipolar 
component of the flux asymmetries can be observed in flgure 
(c): depending on their propagation direction the traces look 
either bright or dark. 
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FIG. 2. Time evolution of heat flux proflle demonstrating 
turbulence movement in responce to a linear shear flow with 
different sign for the lower and upper half of the plots. 



TABLE I. Mean GAM energy production due to Reynolds 
stress and Stringer- Winsor effect, flow intensity (v^) and 
anomalous pressure transport (pVr) for varying parameters. 
The units are the turbulence units described in reference H. 
In these units, the diamagnetic velocity is equal to ad- The 
parameters are consistent with the physical reference param- 
eters given in the text, except for an altered temperature 
Ti — Te — T. For Ud = the ions have been assumed cold. In 
the line marked with an * the Stringer- Winsor drive has been 
switched off, eliminating its braking force in the ballooning 
regime. 

«d Vi T (eV) {pvr) {v%) Reynolds drive Stringer drive 
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